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ABSTRACT 



We compare our calculations of the frequencies of the fundamental g, c, and 
p-modes of relativistic thin accretion disks with recent observations of high fre- 
quency QPOs in X-ray binaries with black hole candidates. These classes of 
modes encompass all adiabatic perturbations of such disks. The frequencies 
^ \ of these modes depend mainly on only the mass and angular momentum of 

the black hole; their weak dependence on disk luminosity is also explicitly in- 
dicated. Identifying the recently discovered relatively stable QPO pairs with 
the fundamental g and c modes provides a determination of the mass and an- 



H ; gular momentum of the black hole. For GRO J1655-40, M = 5.9 ± I.OM0, 

J = (0.917 ± 0.024)GM^/c, in agreement with spectroscopic mass determina- 
tions. For GRS 1915+105, M = 42.4 ± T.OMq, J = (0.926 ± 0.020)G'MVc or 
(less favored) M = 18.2 ± 3. IM©, J = (0.701 ±0.043)GMVc. We briefly address 
the issues of the amplitude, frequency width, and energy dependence of these 
QPOs. 

Subject headings: accretion, accretion disks — black hole physics — gravitation 
— relativity — X-rays: binaries 
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Introduction 



A few years ago, we predicted that the angular momentum of the candidate black hole in 
the X-ray binary GRO J1655-40 is approximately 93% of maximum (Wagoner 1999). (The 
dimensionless angular momentum parameter a = cJ/GM'^ is less than unity in absolute 
value.) This was based upon the identification of the 'stable' 300 Hz feature with the 
fundamental g-mode in its accretion disk (Perez, Silbergleit, Wagoner & Lehr 1997), plus 
the mass determination of 7.02 ± 0.22 Mq from observations of its companion by Orosz & 
Bailyn (1997). As we shall see below, the fundamental c-mode should then be at a frequency 
near 450 Hz, in agreement with the recent discovery of another apparently stable feature in 
its power spectrum by Strohmayer (2001a). A second relatively stable high frequency QPO 
has also been recently seen in another galactic microquasar, GRS 1915+105 (Strohmayer 
2001b). 

For the past ten years, our group has investigated the normal modes of oscillation of 
standard (geometrically thin) equilibrium models of black hole accretion disks. The early 
work of various groups has been reviewed by Kato, Fukue & Mineshige (1998) and Wagoner 
(1999). We perturb the optically thick fully relativistic models of Novikov & Thorne (1973); 
Page & Thorne (1974). We have usually considered accretion disks which are barotropic 
[p = P(p)]) producing a vanishing buoyancy frequency. The results below are relatively 
insensitive to the value of the viscosity parameter, here taken to be a = 0.1. 

The key frequencies, associated with free-particle orbits in the disk, are the rotational 
[f2(a, r)], vertical epicyclic [r2^(a, r)], and radial epicychc [/t(a, r)] angular frequencies, respec- 
tively. The effective inner edge of the disk is close to the radius ri{a) of the last stable circular 
orbit, where /t(a, Tj) = 0. For a > 0, which we shall henceforth assume, Q{r) > Q±{r) > K,{r). 

We stress that linear combinations of these modes should describe all adiabatic per- 
turbations of such disks. There are essentially three classes of modes, designated g, c, and 
p. We shall focus mainly on the fundamental mode within each class, corresponding to the 
minimum number of nodes in the axial, radial, and vertical components of the eigenfunction 
Sp = A{r, z) exp[i(m0 + at)] of the pressure (and any other) perturbation. The correspond- 
ing mode numbers are m (axial), j (vertical), and n (radial). One would expect that most 
mechanisms of excitation (such as turbulence) would produce the highest amplitudes for the 
lowest modes numbers, where the wavelength of the mode best matches the dominant length 
scale of the (turbulent, etc.) excitation. In addition, the resulting luminosity fluctuations 
would tend to increasingly cancel for higher mode numbers. 

The g (inertial-gravity)-mode oscillations (Perez, Silbergleit, Wagoner & Lehr 1997) can 
be characterized by a restoring force that is dominated by the net gravitational-centrifugal 
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force. The axisymmetric {m = 0) g-modes are centered at the radius Vg where the radial 
epicychc frequency k, achieves its maximum. The eigenfrequencies of the modes with only a 
few nodes in the radial eigenfunction are slightly below k,q = n{rg). 

The c (corrugation)-modes which have so far been investigated (Silbergleit, Wagoner & 
Ortega-Rodriguez 2001) are non-radial {m = ±1) vertically incompressible waves near the 
inner edge of the disk, that precess around the angular momentum of the black hole. They 
do not exist in counter- rotating (— 1 < a < 0) disks. Their fundamental frequency coincides 
with the Lense-Thirring frequency (produced by the dragging of inertial frames generated 
by the angular momentum of the black hole) at the outer boundary of their capture zone, 
in the appropriate slow-rotation limit. 

The p (inertial-pressure)-mode oscillations can be characterized by a restoring force 
that is dominated by pressure gradients, as for ordinary acoustic modes. The fundamental 
p-modes exist near the inner and outer edges of the disk, but higher modes can cover a 
greater area. The frequency of the inner edge fundamental p-mode is comparable to that of 
the fundamental g and c-modes, and much greater than that of the outer edge p-mode, so 
only it will be considered below. Details will be presented by Ortega-Rodriguez, Silbergleit 
& Wagoner (2001). 

Unless otherwise indicated, we take c = 1, and express all distances in units of GM/c^ 
and all frequencies in units of c" /GM (where M is the black hole mass). 



2. Properties of the Fundamental Modes 
2.1. g— modes 

The axisymmetric g-modes are trapped in the region of the disk where the eigenfre- 

1 /2 

quency \(j\ < k. The width Ar of this region [centered on rg{a)\ is proportional to cj ■ 
The speed of sound Cs{r, z) is proportional to the disk luminosity L where radiation pressure 
dominates. This should be valid (for typical luminosities) because r^ is near the radius where 
the disk temperature peaks. For < a < 0.95 and small values of the mode numbers j and 
n, the width 

Ar^3.3(L/L^,,)i/% (1) 

while Tgi^) = 8 and rg(0.95) = 3 [see Perez, Silbergleit, Wagoner & Lehr (1997), Figure 2]. 

The eigenfrequencies of the axisymmetric g-modes with small values of the vertical (j) 
and radial (n) mode numbers are 

ag = {1 - ej„)/to , ejn ~ Ejn{a)L/LEdd ■ (2) 
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The small correction term ej„ depends upon the adiabatic index F as well as the usual 
parameters M, a, the viscosity parameter a, and the luminosity L. For relevant values and 
n ^ 3 <l, e ^ h/rg, giving Ej„(0.0) = 0.077 and ^j„(0.9) = 0.11. Thus the frequencies 
should drop slightly as the luminosity increases (as long as the disk remains thin). In 
particular, 

dlogag/dlogL = -Ejn{a)L/LEdd (3) 

when ejn <^ 1. 



2.2. c-modes 

The fundamental {\m\ = j = 1) c-modes are trapped within Vi < r < r^. For radial 
mode number n = 0, its width can be approximated by the formula 

r,-r, = iroa-^ni -«)'"% (4) 

where K^ = 0.093, Ki = 0.79, and K2 = 0.20 for M ~ IOM0 and L ~ O.lL^dd- So the 
mode has a significant radial extent only for slowly rotating black holes, although this extent 
increases with increasing radial mode number. 

The eigenfrequencies of the fundamental c-modes are also relatively insensitive to the 
properties of the accretion disk (except for very slowly rotating black holes). They are also 
related to the basic frequencies mentioned above, in this case by 

ac = ^{rc) - 

We have found that 

d log ac/d log L 

with ir(0.003) = 0.51 and K{0.1) = 0.044 for M 
decrease slowly with increasing values of n. 



2.3. p— modes 

The fundamental {m = j = 0) p-modes that can produce the high frequencies of interest 
here also exist near the inner edge of the disk, Vi < r < Vp, where \a\ > n. We take the 
modified speed of sound 

cl{r) = Cs{r, z = 0)/[U'gl!.'] oc (r - r,)^ (7) 



fi±(re) . 


(5) 


= -K{a) , 


(6) 


M ~ IOMq and n = 


= 0. The frequencies 
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for r close to Vi, where U" is the four- velocity of the fluid and grr is a metric component. For 
the standard disk model, n = 2/5. However, magnetic stresses in the gas that is spiraling 
toward the black hole (Hawley & Krolik 2001) may torque the inner edge of the disk enough 
to produce /x = 0. The fraction of the specific angular momentum at r, that is eventually 
captured by the black hole is usually denoted by p. [Thus /i(/9 = 1) = 2/5 and /i(/9 < 1) = 0.] 
For < fx < 1/2, we find that 

r,-r, = [CMq-'c:{r,)f\ (8) 

where q'^{a) = (dK'^/dr)^ and C„ ~ {2n + l)[2/(r + 1)]"*^/^ for any value of /i. Since typically 
Tp — Ti<^ 1, we expect gas pressure to dominate, giving F = 5/3. For a < 0.9, this equation 
can approximated by 

rp-r, = D„(/?,L)K(AL)-a]. (9) 

For instance, DQ{l,Q.lLEdd) = 0.041 and al{l,0.1LEdd) = 1-22, while Do{0.95,0.1LEdd) = 
0.24 and ao(0.95, O.lLgdd) = 1.18. This trapping region size increases (roughly as n^^^) for 
increasing values of the radial mode number n (like the c-modes). 

The eigenfrequencies of the fundamental inner p-modes have a similar dependence on 
the speed of sound. In terms of the quantities defined above, the basic relation a"^ = n^ij-p) = 
q^{rp -Ti) gives 



a. 



p 



[Cn{fi)q'c:{r,)Y/' , (10) 

so (r-p — ri)ap = Cn{lj)c*g{rp). This frequency increases (roughly as n^^"^) for increasing values 
of n. Since Cg oc [L"^ /aMY^^^ near the inner edge, 

dlogap/c^logL = 1/15. (11) 



3. Comparison with Observations 

In Table 1 we indicate the radial location and extent of each mode, for two values of the 
angular momentum parameter a. In Figure 1 we plot the fundamental eigenfrequencies of 
each mode as a function of a. All the mode frequencies are essentially proportional to 1/M. 

Since all of the mode frequencies are relatively weak functions of the disk luminosity 
(and other disk properties such as the viscosity parameter), we shall focus our attention on 
the two black hole candidates which exhibit more than one 'stable' high frequency quasi- 
periodic oscillation (QPO) in their power spectrum. They are the microquasars GRO J1655- 
40 (Remillard et al. 1999b; Strohmayer 2001a) and GRS 1915+105 (Morgan, Remillard & 
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Greiner 1997; Strohmayer 2001b). In Table 2 we present their QPO frequencies /, qual- 
ity factors Q = f /Af, and fractional amplitudes. The amplitudes are usually increasing 
functions of photon energy. 

We can obtain numerical estimates of the dependence of the frequency / on the observed 
energy flux JF from equations (3), (6), and (11). For the p-modes, (ilog//(ilogjF = 0.067. 
Note that the frequencies of the g and c-modes decrease with increasing flux. We see that 
\d\ogf/d\ogJ-'\ < 0.1 for them ii L < LecLcl (necessary for our thin disk assumption) and 
a > 0.03. Only the c-mode can produce a greater variation, for very slowly rotating black 
holes. Unfortunately, there is little information available on the disk (thermal) flux when 
each QPO was observed. So all we can do is use these limits to eliminate from consideration 
any QPOs which violate them. For instance, if we assume that the thermal flux JF varies by 
a factor of < 5, then we require that the frequency / vary by < 16%. This allows the shift 
in frequency of about 4% observed in GRS 1915+105, if the flux was larger when the lower 
frequency was observed. 

From Table 1, we see that the p-mode probably occupies a smaller region than the other 
modes. In addition, its oscillations have very little emergent flux to modulate since they are 
located near the inner edge of the disk. We will therefore focus on the g and c-modes. Then 
the pair of frequencies observed in each source allow (at most) two possible choices for M 
and a (with a spread from unknown properties of the disk), as seen from Figure 1. 

For GRO J1655-40, choosing the lower frequency (295 Hz) to be that of the g-mode [in 
keeping with our earlier choice (Wagoner 1999)] gives M = 5.9 ± l.OM©, a = 0.917 ± 0.024. 
This predicted mass is consistent with the value M = 6.3 ± O.5M0 obtained from spectral 
observations of the companion star (Shahbaz et al. 1999; Greene, Bailyn & Orosz 2001). 
(From Figure 1, but for 0.1 < L/LEdd ^ 1-0, the corresponding p-mode frequency is then 
118—193 Hz for f3 = 0.95 and 50.4—87.6 Hz for (3 = 1.00.) Identifying instead the lower 
frequency with that of the c-mode gives an unacceptably smaller mass (~ 2.5M0). 

However, there is an apparent problem with this assignment of the modes. Since the 
g-mode occupies a much larger region of the disk (Table 1, a = 0.9), one might expect 
that the resulting amplitude of the luminosity modulation would be greater than that of 
the c-mode if they are excited in the same manner. This disagrees with the amplitudes in 
Table 2. However, since the c-mode is virtually incompressible, it should mainly modulate 
the coronal photons via reflection from its changing projected area; whereas the g-mode 
modulates the internal properties of the disk (Nowak & Wagoner 1993; Nowak, Wagoner, 
Begelman & Lehr 1997). This matches the fact that the higher frequency QPO has the 
harder energy dependence (Strohmayer 2001a), since the coronal (power law) photons are 
more energetic than those emitted directly by the disk. In addition, the viscous damping 
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(or growth) rate of the c-mode should be less than the viscous growth rate of the g-mode 
(Ortega- Rodriguez & Wagoner 2000), producing a larger Q value (which agrees with Table 
2). Of course, nonlinear it ies may play an important role. The excitation of higher mode 
numbers, as well as radiation leakage, can also contribute to the frequency width (Nowak, 
Wagoner, Begelman & Lehr 1997). 

For GRS 1915+105, our earlier choice (Nowak, Wagoner, Begelman & Lehr 1997) of the 
g-mode for production of the 67 Hz QPO gives M = 18.2 ± 3.1Mq, a = 0.701 ± 0.043. (The 
corresponding p-mode frequency is then 34.6 — 56.8 Hz for f3 = 0.95 and 14.5 — 25.2 Hz for 
f3 = 1.00.) If it is produced by the c-mode, M = 42.4 ± 7.OM0, a = 0.926 ± 0.020. (The 
corresponding p-mode frequency is then 16.3 — 26.5 Hz for (3 = 0.95 and 7.00 — 12.1 Hz for 
(3 = 1.00.) This second choice is supported by the claim of Strohmayer (2001b) that the 
dependence of the properties of QPOs on their frequency is in the same direction in both 
sources, as moderately indicated by the data on the quality factors and amplitudes of the 
QPOs in Table 2. 

Not unexpectedly, the observed power spectra of both sources show no 'stable' features 
at the p-mode frequencies indicated above. An example of another mode of interest is the 
m = 1 g-mode, whose frequency is about a factor of 4 greater than that of the m = 
(fundamental) g-mode (Perez, Silbergleit, Wagoner & Lehr 1997). The large inclination 
angle of about 70 degrees for both of these microquasars could make this nonaxisymmetric 
mode (like the c-mode) visible. However, existing power spectra are often too noisy to be 
able to detect such higher frequencies. 

There is not yet evidence that any of the high-frequency QPOs in other black hole can- 
didates [XTE J1550-564 (Remillard et al. 1999a; Homan et al. 2001; Miller et al. 2001), XTE 
J 1859+226 (Cui, Shrader, Haswell & Hynes 2000), and 4U 1630-47 (Remillard & Morgan 
1999)] are stable to the extent required if they are produced by these modes. Some are known 
to vary significantly, while others have only been detected during a single observation. 

This work was supported by NASA grant NAS 8-39225 to Gravity Probe B. 
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Table 1. Mode Trapping Regions* 





P 


Radial location 
a = 0.1 a = 0.9 


Radial extent 


Mode 


a = 0.1 


a = 0.9 


9 


— 


7.6 


3.4 


1.0 


1.0 


c 


1.00 


5.7 


2.3 


0.56 


0.064 


P 


0.95 


5.7 


2.3 


0.26 


0.067 


P 


1.00 


5.7 


2.3 


0.046 


0.013 



*For L ~ O.lLEdd, a ~ 0.1, and M ~ IOMq. Units 
are GM/c'^. 



Table 2. QPO Properties* 



Source 



Frequency (Hz) Q Amplitude (%) 



GRO J1655-40 




295 ±9 


^ 4 


0.81 ± 0.11 


GRO J1655-40 




449 ± 5 


~ 23 


4.8 ± 0.6 


GRS 1915+105, 


AOl 


66.4 ± 1.5 


~ 20 


6 ± 1 


GRS 1915+105, 


A02 


69.2 ± 0.2 


~ 20 


1.9 


GRS 1915+105, 


A02 


41.5 ± 0.4 


r^ 8 


2.4 



*Froni Remillard et al. (1999b), Strolimayer (2001a), Morgan, 
Remillard & Greiner (1997), and Strolimayer (2001b). AOl and 
A02 refer to different epochs of observation. 
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Fig. 1. — The dependence of the fundamental frequencies of the three diskoseismic modes 
on the angular momentum of the black hole. The spread of each corresponds to the range 
0.01 < L/LEdd ^ 1-0, and is relatively insensitive to the choices M/Mq ~ 10 and a ~ 0.1. 
The upper band for the p-mode corresponds to /i = with (3 = 0.95, and the lower band to 
/i = 2/5 (no torque, P = 1.00). 



